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ABSTRACT 

This work reports on the apphcation of the Eulcrian perturbation theory to a re- 
cently proposed model of cosmological structure formation by gravitational instability 
(Dommguez 2000). Its physical meaning is discussed in detail and put in perspective 
of previous works. The model incorporates in a systematic fashion corrections to the 
popular dust model due to multistreaming and, more generally, the small-scale, viri- 
ahzed degrees of freedom. It features a time-dependent length scale L{t) estimated to 
be L/tq ^ i^oit) is the nonlinear scale, at which (5^) = 1). The model provides a 
new angle on the dust model and allows to overcome some of its limitations. Thus, the 
scale L{t) works as a physically meaningful short-distance cutoff for the divergences 
appearing in the perturbation expansion of the dust model when there is too much ini- 
tial power on small scales. The model also incorporates the generation of vorticity by 
tidal forces; according to the perturbational result, the filtered vorticity for standard 
CDM initial conditions should be significant today only at scales below ^ 1 h^^Mpc. 
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1 INTRODUCTION 

The evolution of cosmological structures by gravitational in- 
stability is usually modelled with the popular dust model 
(Peebles 1980; Sahni & Coles 1995): it is simply the hydro- 
dynamic Eqs. for a fluid under the influence of no other force 
but its own gravity (no pressure, no viscosity, no heat flows). 
The usual justiflcation of this model is that (long-range) 
gravity is the overwhelmingly dominant force on the cosmo- 
logical large scales and that, while the gravitational instabil- 
ity is not too advanced, the matter distribution looks like a 
continuum. It is found, however, that the success of the dust 
model extends beyond the expectations raised by this ar- 
gumentation, as evidenced by the comparison with N-body 
simulations of the perturbation theory predictions, both Eu- 
lerian (Juszkiewicz et al. 1993; Bernardeau 1994; Gaztanaga 
& Baugh 1995; Bernardeau 1996) and Lagrangian (Buchert 
et al. 1994; Melott et al. 1995; Bouchet et al. 1995; Weiss 
et al. 1996): some perturbation predictions hold even close 
to the nonlinear scale ro, when the matter distribution does 
not look homogeneous any more. But the dust model has of 
course its own shortcomings, notably that the solutions typ- 
ically develop singularities: this already happens in the first 
order Lagrangian solution (Zel'dovich 1970; Buchert 1989; 
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Moutarde et al. 1991) (the singularities are generically sheet- 
like — 'pancakes' in the cosmological literature); the absence 
of any dissipative term in the dust model Eqs. suggests that 
this is a property of the model and not an artifact of the 
perturbation expansion. Nevertheless, this problem could be 
satisfyingly (when compared to simulations) tackled by a 
phenomenological correction, the adhesion model (Gurba- 
tov et al. 1989) , which adds a viscosity to the dust model. 

Recently, I proposed a novel approach (Dommguez 
2000; Paper I hereafter). The new model is characterized by 
a correction term to the dust model featuring a length scale 
L. It relies on a formal expansion in powers of L which I call 
the small-size expansion (SSE)^. I showed in Paper I that the 
correction behaves like an effective viscosity for the generic 
sheet-like collapse configuration; in this sense, it represents 
a derivation of a generalized adhesion-like model which con- 
tinues the line of research of (Buchert & Dommguez 1998; 
Buchert et al. 1999). In that work I also raised two ques- 
tions, which will be the goal of this paper: a better physical 



t This anticipates the discussion in Sec. 5; in Paper I the denom- 
ination was the 'large-scale expansion', because it can be viewed 
also as a formal expansion in powers of k in Fourier space. But 
this designation is unfortunate, since it can be mistaken for the 
large-scale cosmological expansion. 
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understanding of the length L and the possibihty that the 
corrections may act as a source of vorticity. 

The analytical approach is hindered by the nonlinear 
nature of the dust model and the corrections: in Sec. 2 I 
apply the Eulerian perturbation expansion in the nonlin- 
earities, which will allow a direct comparison to the results 
for the dust model obtained with the same technique. The 
calculations are mathematically a bit more involved, but 
there are no essential differences. The skewness of the den- 
sity contrast is computed to lowest order in Sec. 3, since this 
example provides a test case for comparison of the SSE pre- 
dictions with the dust-model ones. Sec. 4 studies the SSE 
prediction to lowest perturbational order for the vorticity of 
the peculiar velocity. This is most interesting because the 
vorticity vanishes in the dust model and is thus entirely due 
to the correction. Sec. 5 provides a detailed discussion of the 
physical meaning and interpretation of the SSE and the scale 
L. Finally, Sec. 6 summarizes the conclusions and points out 
several possible lines of future research. App. A briefly re- 
peats the derivation of the SSE in Paper I. App. B collects 
mathematical manipulations. 



2 CORRECTIONS TO DUST AND 
PERTURBATION THEORY 

The model with corrections to dust which I proposed in Pa- 
per I introduces a comoving length scale L through smooth- 
ing of the microscopic mass and peculiar-momentum density 
fields. Then it is assumed that the coupling of the small-scale 
(< L) degrees of freedom with the large-scale (> L) ones is 
weak. This allows a certain expansion in powers of L (this 
is the SSE) to be truncated. App. A contains a brief deriva- 
tion of the model, slightly generalized over the one in Paper I 
by allowing for a time-dependent smoothing length L{t). It 
yields the following Eqs. in standard comoving coordinates 
for the mass density field q{x, t) and the peculiar velocity 
field u{x,t) to first order in the SSE (5, = d/dxi, and a 
summation is implied by the repeated index i): 

at a L 

= -AHqu + OW--V ■ (ouu) + 
at a 

+ BLU{Vg ■ V)w - iv • [Qidiu)id^u)] + yV^gu] 



V • w = —4:nGa{g — gt), 

V X to = 0, 



(1) 



where i? is a constant of order unity fixed by the shape 
of the smoothing window. When L is formally set to zero, 
Eqs. (1) reproduce the dust model and this means physically 
that the structure below the scale L is dynamically unimpor- 
tant. The terms proportional to L have a simple geometrical 
meaning, as discussed in App. A. The other two correction 
terms represent the influence of the small-scale structure 
on the dynamical evolution of the large scales: (V^i ■ V)'w 
is a tidal correction to the 'macroscopic' gravitational field 
to; V • [g{diu){diu)\ is a correction to the advective term 
V • {guu) due to the velocity dispersion induced by shears in 



the 'macroscopic' velocity field u. Conceptually, these cor- 
rections have the same origin as pressure, viscosity, heat 
conduction in usual hydrodynamics, namely, the dynamical 
effect of the microscopic degrees of freedom neglected in the 
coarsed description. But of course, their form is very differ- 
ent because of the qualitatively disparate underlying physics 
(long vs. short range dominant interaction). 

The nonlinear character of Eqs. (1) makes it difficult to 
learn about the influence of the corrections without resort to 
some approximation. In this work I will employ the Eulerian 
perturbation expansion: introduce a bookkeeping parameter 
e by writing 



(2) 



where 5 :— {g/ gb) — ^ is the density contrast. The parameter 
e is recognized to count the degree of nonlinearity. Pertur- 
bation theory consists in looking for a formal solution as an 
expansion in powers of e which is truncated at some order; 
at the end of the calculations, one sets e = 1 to recover the 
original problem. This has been a widely used tool to study 
the dust model (Bernardeau et al. 2002). 

This procedure assumes that the departures 5, u from 
the homogeneous expanding background are small. In the hi- 
erarchical scenario I will consider, this is satisfied when the 
fields are observed on sufficiently large scales. This implies in 
turn that the perturbational solutions should be employed 
in principle only to estimate quantities which axe already de- 
fined by smoothening over some sufficiently large scale. But 
this is not enough, because the nonlinearities couple widely 
different scales. Consequently, some of these perturbation- 
ally computed quantities may exhibit a short-distance, ul- 
traviolet (UV) divergence within the dust model [L — »• 0); 
this divergence is however naturally regularized in the con- 
text of the SSE by the finite scale L. This will be seen in 
next Sees, with examples. 

The perturbation expansions (2) axe inserted in 
Eqs. (1). The initial conditions at time tin read 

5l{x,tin) = 5{x,tin), 5x{x,tin)=Q, A > 2, 

and similarly for u. To lowest order the linearized version of 
Eqs. (1) is found: 

^ = -iv-Mi-FSLiV^5i, 
at a 

^ = -Hui -h toi -h BLLV'^ui, 
at 

V • = —A^GagbSi, 

V X wi = 0. 



(3) 



These Eqs. differ from the well-known linear dust model only 
in the terms oc L. The solutions axe more easily written down 
for the Fourier transform of the fields (to avoid burdening 

the notation, the same symbol will be used for a quantity 
and its Fourier transform. Which one is meant should be 
clear from the argument of the function and the context). 



dx e 



{x). 



In the long-time limit (t ^ tin, but still within the validity 
regime of the perturbation expansion), the solution reads 
{Lt = L{t),k=\k\y. 
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Si{k,t)=bte Smic{k,tin), 



(4) 



ui{k,t) = -atbt^e Smic{k,ti„), (5) 

where bt is the growing mode of & + 2Hb — ATvGgbb = 0, with 
the normahzation a{tin) = b{tin) = 1, and 

^mic{k^tin^ = 6^ ' S(k.tin^ 

is the unsmoothed initial density contrast, Eqs. (A3). (This 
identification is rigourously true only for a Gaussian window; 
otherwise, it is part of the approximation behind the SSE, 
see end paragraph in App. A). Compared to the usual lin- 
ear solutions, Eqs. (4-5) exhibit the explicit time-dependent 
smoothing. This spoils the separability of the time and 
wavevector dependences and renders the calculations a bit 
more messy. 

Knowing the solutions (4-5), the first-order nonlinear 
corrections obey a set of linear, inhomogeneous, partial dif- 
ferential Eqs.: 

^ = - iv • «2 + BLLV''62 - -V ■ (JiMi), 
at a a 

= —Hu2 + W2 + BLL'V^U2 — — (tti • V)mi + 
ot a 

+ BL^ |(V5i • V)w;i - iv • [{diui){diui)] + 



-2-(VJi • V)ui 

1^ 



(6) 



V • k;2 = —4:7vGaQbS2, 

V XW2 = 0. 

These Eqs. will be studied in the following Sees. As re- 
marked, the time-dependence of L complicates a little the 
calculations. However, the qualitative features of the results 
to follow axe the same as if L = 0. 



3 SKEWNESS OF THE DENSITY CONTRAST 

As an example of how to proceed, Eqs. (6) are solved for 
$2 in this Sect, and the skewness of the density contrast is 
computed to leading-order. The terms oc L are eliminated 
by introducing auxiliary variables: 



A2{k,t) := <52(fe,t)e3^^*'= , U2ik,t) := U2ik,t)e^ 
Eqs. (6) then read: 

dA2 1., rr 

—— = -ik -1/2+ Qa, 
at a 

dU2 rrrr . ^ A ^ 

-Q^ = -HU2 - A-KGaQbj^ A2 + Qu, 
with the sources: 



(7) 



QA{k,t) := 



e2 * 
at 



dq 



(2^) 



k ■ ui{q,t)]6i{k — q,t), 



Quik,t) :=e5«^?'='y' ^i^l.[iq.mik-q,t)]miq,t) + 



+ BLtlq ■ (k - q)] 



to 

4'KGatQbSi{q,t)6i{k - q,t)-^ - 



1 2L 

[ik ■ ui{q,t)]ui{k — q,t) ^Si{k — q,t)ui(q,t) 

at J-it 

Eqs. (7) now look the same as those of the dust model except 
for the especific form of the sources. The solution for A2 with 
the initial condition A2{k,tin) = is 



A2{k,t) = bt f dr Tdr' 



Q{k,t) := 2HtQA{k,t)+QA{k,t) 

dq 



Q{k,r') 



-ik ■ Qjj{k,t) = 



k • q 



at 

S ■ (a f- )S ■ (k-a f- \ABLUk^-q^-\k-q\^) 



[1 - BLl q ■ (fe 



'inGQb - + 



k-{k-q) 



where the linear solutions (4-5) have been inserted. To sim- 
plify the solution further, we assume an Einstein-de Sitter 
background, bt = at = {t/Un)^^^- Then, one can write fi- 
nally: 



52{k,t) 



= [ ^ 
J (27r: 



dq_ 



X Gsik - q,q)J^s{k - q,q,t), 

with the kernel 

q-{k + q) 



(8) 



G6ik,q) := 



3 k-{k + q) 
2 fe2 



and the function 

T,{k, q, t) := a*e-5^^?l'=+-l' J^' dr £ dr' (^^^ ^ x 

x[l-BLl,q.k] e^BLl,(\k+q\--e-k-)_ 

The corrections to dust are collected in the function Ts- In 
the dust limit, L ^ 0, it becomes a function of time only: 
(2/7)a? (in the long-time limit t S> tin)- 

The second-order solution can be employed to esti- 
mate the skewness of the density contrast, defined as S3 := 
where (■ • •) denotes an ensemble average over ini- 
tial conditions. When the initial density contrast is Gaussian 
distributed, (Sf) = by Eq. (4). Hence, the skewness of the 
evolved density contrast is generated solely by the nonlinear 
evolution of the gravitational instability and this motivates 
the interest in its computation. 

The variance and the skewness of the density contrast 
smoothed on a scale R read: 



a\R,t) = J 
S3{R,t) 



^{S{k,t)5{q,t))W{Rk)W{Rq), 

1 



(9) 



(T^{R,t) J (27r) 

xW{Rk)W{Rq)W{Rp), 

where W{-) is the Fourier transform of the smoothing win- 
dow. As mentioned previously, the use of the perturbation 
theory solutions (2) requires R to be large enough. Further- 
more, initial Gaussian inhomogeneities are assumed, char- 
acterized by the power spectrum 



{5mic{k,tin)5mic{q,tin)) = {2nfP{k)S^^\k + q) 



(10) 
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With the perturbation solutions (4) and (8), the variance to 
leading order reduces to the linear solution, 

dfe 



ae (R, t) = at 



(2^) 

while the skewness reads 



^P{k)\W{Rk)\e 



2 -BL^k^ 



S3{R,t) = 



af{R,t) J (2n 



dk dq 



P{k)P(q)g,{k,q) X 



(11) 



(12) 



X J^s{k, q, t)WiRk)WiRq)W{R\k - q\). 

Comparison with N-body simulations has shown that the 
dust prediction {L = 0) is very good (Juszkiewicz et al. 
1993), and not only for very large R but also oven close to 
the nonlinear scale ro, defined by the condition a{ro, t) = 1, 
where perturbation results would be expected to break 
down. In fact, the limit L — > can be taken inside the inte- 
grals in Eqs. (11) and (12), since the smoothing on scale R 
already assures UV convergence. Taylor-expansion provides 
a correction to the dust prediction of order (L/R)^, and this 
suggest that the scale L in the SSE is smaller than ro. 

The correction to dust is barely noticeable in the 
leading-order skewness because the scale R already works 
as a short-distance cutoff. But this is of course not the case 
at higher perturbational orders or with other measurable 
quantities. The simplest example is (t^{R): it can be argued 
(Valageas 2002) that the perturbational correction in the 
dust model always exhibits an UV divergence at a suffi- 
ciently large order, in spite of the smoothing over a scale 
R well in the linear regime. Indeed, the next-to-leading per- 
turbational contribution to <y^{R) (of order e'^, which re- 
quires going to third order in the perturbation expansion) 
exhibits the following scaling when P{k) (x fc" (self-similar 
power spectrum) within the dust model (Lokas et al. 1996; 
Scoccimarro & Frieman 1996b): 



\R,t} 



2 

-ere 









R 











-3 <n< -1, 



2{n+3) 



-(n+1) 



-1 < n, 



with an UV cutoff L <^ R. If n < —1, the smoothening over 
the large scale R is enough to avoid short-distance diver- 
gences and the scale-dependence can be guessed from self- 
similarity arguments. However, if n > —1, there is too much 
small-scale power and the non-linear couplings require an 
additional UV cutoff L. 

In the next Sec, an example is studied where the de- 
viations from the dust model are particularly evident: the 
vorticity of the velocity field. 



4 VORTICITY OF THE PECULIAR 
VELOCITY 

The dust model Eqs. lack a source of vorticity, u; := V x u: 
it can arise only from an initial vorticity, which is neverthe- 
less strongly damped in the linear regime. This tiny initial 
vorticity may be amplified in high-density, collapsing regions 
(Buchcrt 1992; Barrow & Saich 1993), but then caustics de- 
velop (multistream flow, the field u becomes multivalued) 
and the dust model itself breaks down. Thus, the correc- 
tions to the dust model are most relevant in this case. As 



I will show, Eqs. (1) do indeed generate vorticity, and the 
leading-order perturbational contribution will be computed. 

Let us take the curl of Eq. (6) for U2 and use that 
wi = by virtue of the solution (5), to get 

dui2 
~dt 



+ Hu2 - BLL\/U2 = 



(14) 



BL^V X { {V5i ■ V)wi - -V • l{diui){diui)] + 



+ 2-(V5i • V)wi 



(15) 



Upon introducing the traceless shear tensor of the peculiar 

velocity field, 

~ ^{diUj + djUi) - i(V ■ u)Sij, 

and the tidal tensor of the peculiar gravitational accelera- 
tion, 

Eij ■- ^{dtWj + djWi) - i(V • w)5ij, 

one can write the source term alternatively as 



BL'V X El • V5i 



-V- 



^(V • u\_)a\ + ai • (Ti 



-h2-(V5i ■ V)ui 

Since ai is in fact proportional to Ei by Eq. (5), the ultimate 
dynamical source of vorticity is the gravitational tidal force, 
whose effect is manifested directly (first term) and indirectly 
through the induced velocity dispersion (second and third 
terms). As already remarked, the fourth term has a purely 
geometrical origin. 

Inserting the linear solutions (4-5) in Eq. (15) and 
Fourier transforming, this source reads 

(13) Q^{k,t) = BL"^^ I ■jP-.Smiciq,tin)Smic{k - q,tin) X 

at J (2'!vy 



(fc ■ qjjk X q) ^_ 1 BL2(^2_^n,_g|2) 



(16) 



JF„(t) := {abf 



1 -I- AnGgt 



+ 2 



bL 
bL 



Eq. (14) can be integrated immediately for the Fourier trans- 
form of the vorticity with the initial condition u!2{k,tin) = 
0: 



u;2(fc,t) = Tdr ^e-^«(^*-^?"=^QJfe,r). 



(17) 



An estimate of the strength of the vorticity on a scale 
R is provided by its variance, 

dfe dq 



at{R,t) 



-{w{k,t) ■ w{q,t)}W{Rk)W{Rq). 



(27r)6 

Inserting the perturbation solution (16-17) and simplifying: 
dfe dq 



'yt{R,t) 



'-P{k)P{q)gUk,q)\W{R\k-q\ 



i2^r 

x|H(fe.9,t)|^e-^^?"=-''l' 
with the positive, symmetrized kernel 



(18) 
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1 



and the function 

1 /■* 2 



(19) 



To gain insight into this result, it will be worked out for 
the particular case of a self-similar model: an Einstein dc 
Sitter background and an initial power spectrum P{k) (x fc" 
are assumed. It remains to specify the smoothing length 
L{t). In a self-similar model, the only relevant scale is ro{t) 
and the natural choice is L{t) — Aro(t), where A is a time- 
independent proportionality constant. This choice will be 
supported by the discussion of the physical meaning of the 
SSE in the next Sec. 

The details of the mathematical analysis are collected 
in App. B. The physically interesting results can be summa- 
rized in the asymptotic behavior of a^i in the limit Lt 
and t:$>ti„: 

t-3 









2 








.R. 




R 


d{t) 






2 








R. 




R 



-3 < n < - 



2' 



(20) 



n+3 



3 

-2<n, 



up to a dimensionless factor which depends only on the spec- 
tral index n. The lower bound — 3 < n prevents an infrared 
(IR) divergence. The crossover at n = —3/2 is related to 
the UV behavior of the integrals, and we arc in a simi- 
lar situation as with Eq. (13). If n < —3/2, R is enough 
to regularize possible UV divergences, the scale L only ap- 
pears in the prefactor BL^ of the source (15) and the scale- 
dependence of can be guessed by simple power-counting. 
If n > —3/2, however, L is relevant also as UV cutoff. 
Fig. 1 shows the numerical computation of ctw with Eq. (18), 
confirming the asymptotic scaling (20); in particular, that 
the i?-dependence is insensitive to the spectral index when 
n > —3/2. It also shows how well the asymptotic scaling 
is already obeyed when R/Lt ~ 10. In the opposite limit 

— > (known in the literature as the 'unsmoothed' limit), 
(Ta> approaches a nonvanishing constant. 

To address the relevance of the generated vorticity, it 
will be compared to the dust prediction for the divergence 
of the velocity field, ^ := V • u. The variance erg := {6^) on a 
scale R is defined in analogy with Eq. (9). The variance at, 
of perturbational order e*, must be compared with the next- 
to- leading order in a1, denoted by ag. This term scales like 

times Eq. (13); more generally, each term in the perturba- 
tional expansion of ag is proportional to the corresponding 
one in the perturbational expansion of {aaY (Scoccimarro 
& Prieman 1996a; Fosalba & Gaztanaga 1998). Thus 



-3 < n< -- 



5-9 









R 




Lt' 




R 




Lt 




R 



-- <n< -1, 



(21) 



-1 < n. 



up to a factor of order unity. The exponent of the ratio Lt /R 




0.5 1 1.5 
log(ii/Lt) 

Figure 1. Log- log plot of the numerically computed auj as a 
function of R/Lt in the long-time limit for an initial P{k) oc fc" 
and for a Gaussian smoothing window. The lines have the slope 
predicted by the asymptotic scaling (20). 



is in all cases positive. At fixed resolution R, the vorticity 
becomes more and more relevant in time; at fixed time, it 
becomes subdominant on large scales. 

A more realistic CDM initial power spectrum at the 
scales of physical interest can be approximated well by a 
power-law P(k) oc fc" with n w —1, —2. This initial condi- 
tion is set at a redshift z ~ 10^, and so the present epoch 
corresponds to at ~ 10'', well in the long-time regime (sec 
App. B). Realistic deviations from an exact Eiristein-de Sit- 
ter background are expected to aS'cct only weakly the re- 
sults, as happens with predictions by the dust model. As 
will be argued in Sec. 5, the ratio A = Lt/ro{t) is ~ 0.1. 
With a value ro(todajr) ~ 8 /i^^Mpc, the estimate (21) sug- 
gests that the vorticity becomes significant today on scales 
ii w Lt w 1 h~^Mpc; the physical vorticity u)/a has then a 
value of the order of 1-10 times the Hubble constant today, 
as can be immediately read from Fig. 1. This conclusion does 
not disprove the approximation of potential flow at weakly 
nonlinear, large scales {R » ^o(t)), which is invariably done 
both in theoretical (Sahni & Coles 1995) and observational 
(Dckcl 1994) studios. Both the length scale of relevance and 
the magnitude of the vorticity agree with the result by Pi- 
chon & Bernardeau (1999), who studied the growth of vor- 
ticity with the Zel'dovich aproximation at times just after 
the formation of the first caustic. But the SSE offers a more 
straightforward and systematic method than the one em- 
ployed by these authors. 



5 DISCUSSION 

As follows from the derivation in App. A, the SSE relies 
on a dynamically negligible effect of the substructure of the 
coarsening cells of size « L; the only dynamically relevant 
properties of the cells are the mass, oc g, and the center-of- 
mass velocity, it. On the other hand, one may notice that 
Eqs. (A1-A2) for the evolution of point particles under its 
own gravity can be formally written as the dust model Eqs. 
for the fields Qmic and Umic, Eqs. (A3). Therefore, an in- 



© 2002 RAS, MNRAS 334, 435-443 



6 A. Dommguez 



terpretation of the assumption behind the SSE is that the 
evolution of the scales > L corresponds, in the lowest ap- 
proximation, to that of a set of 'effective point particles', 
plus corrections due to the nonvanishing particle size « L 
and their internal structure. This interpretation motivates 
the appellation 'small-size expansion'. In a bottom-up sce- 
nario, it is natural to identify the effective particles with the 
most recently (and thus largest) formed clusters. In fact, 
Peebles (1974; 1980, Sec. 28) had already shown the mutual 
cancellation of the nonlinear couplings to the small scales 
in the Eq. for S when these small scales consist of virialized 
clusters. 

It must be noticed that this provides a different phys- 
ical interpretation of the dust model and extends its va- 
lidity. Indeed, the dust model is usually justified start- 
ing from the ideal fluid model: then the pressure term is 
dropped because self-gravity is dominant on scales much 
larger than the Jeans' length. This argumentation applies 
in top-down ('pancake') models because, at a given time, 
matter is smoothly distributed on small scales and the free- 
fly of particles gives rise to a kinetic pressure: the ideal fluid 
model is a good approximation. This does not hold, however, 
in the opposite limit of a bottom- up ('hierarchical') evolu- 
tion, for the small scales are highly structured. Nevertheless, 
the SSE axgues how the dust model can still be a good de- 
scription when the evolution is observed through the right 
spatial resolution ~ L~^: then particles appear trapped in 
the virialized clusters and to a first approximation, there 
is no free-fiy, no kinetic pressure. In both extreme scenar- 
ios is the dust model the lowest-order approximation, the 
differences show up only in the corrections. 

The Zel'dovich approximation (Zol'dovich 1970) was 
originally devised to address structure formation in top- 
down scenarios (HDM models). It is the lowest-order term 
in a Lagrangian perturbation expansion of the dust model 
(Buchert 1989; Moutardc ct al. 1991). Thus, it came out 
as unexpected and worth remarking that it provided also 
a good description of the large-scale structure in bottom- 
up simulations, particularly if the initial power spectrum 
was first smoothed on a scale ~ ro: this is the truncated 
Zel'dovich approximation (P. Coles ot al. 1993; Mclott ct 
al. 1994; Pauls & Melott 1995). The SSE offers a physical 
explanation for this convergence and, as a novel aspect, im- 
plements in a systematic way the role of the 'truncation' 
(smoothing) and the 'residual' coupling to the neglected 
small scales. 

We see that smoothing is an essential component of the 
dust model itself in a hierarchical scenario. The use of an 
UV-cutoff in the perturbation expansion is not new, but it 
was mainly introduced as a matter of practical purposes; 
e.g., to regularize the divergences in the perturbation ex- 
pansion of the dust model as the initial small-scale power 
is increased (Sec. 3). These divergences are regarded un- 
physical, since N-body simulations show nothing special in 
the observables. This is usually interpreted as a failure of 
the perturbation expansion, assuming the results from sim- 
ulations to be qualitatively right, and so it is conjectured 
that divergences could be cured by resumming the expan- 
sion (Bernardeau 1996). The SSE can be viewed as such a 
resummation: it takes account of a highly nonlinear feature 
— the virialized clusters, and its dynamical effect on the 
evolution of the large scales. 



According to the physical interpretation of the SSE, 
the smoothing length L{t) should scale with ro(i), as the 
natural measure of the size of the largest virialized struc- 
tures. Consequently, the self-similar behavior expected when 
P{k) oc fc" on an Einstcin-de Sitter background is not 
spoiled. Such a cutoff was in fact employed by Jain & 
Bertschinger (1996), and criticized in turn by Scoccimarro & 
Prieman (1996b), who rejected its ad hoc time dependence. 
The latter authors suggest that self-similarity breaks down 
in general in perturbation theory, because the UV cutoff is 
just a property of the initial conditions. The SSE approach 
reconciles self-similarity with perturbation theory by realiz- 
ing that the (iynamically generated UV cutoff L(t) oc r(]{t) is 
itself a defining ingredient of the dust model in a hierarchical 
scenario. 

Structures of size L{t) should behave approximately as 
'particles' (clusters); hence, L{t) must be defined through 
a condition which explicitly tests somehow the degree of 
'dynamical isolation' of the structure below a given scale. In 
Sec. 3 it was concluded that L is different from the nonlinear 
scale To: then, cr^ = 1 is not such a tester. Nevertheless, 
a possible suggestion of how the function a{R,t) could be 
useful is the following: Lot us assume that, in the case of 
a scale-invariant initial power spectrum and an Einstein-de 
Sitter background, a{R,t) is given by the stable clustering 
ansatz (Davis & Peebles 1977) in the nonlinear regime {R <C 
ro). Then, the recently formed clusters, and thus the length 
L, may be identified by the scale 7? at which a{R,t) starts 
following the stable clustering prediction. Fig. 3 in (Colombi 
et al. 1996) and Fig. 1 in (Jain 1997) provide the estimate 
L/ro ~ 10~^. N-body simulations seem to support the stable 
clustering hypothesis (Efstathiou et al. 1988; Colombi et al. 
1996; Jain 1997; Valagoas ct al. 2000). But this point is not 
settled yet and there are evidences that it only holds for 
the recently formed clusters (Munshi et al. 1998; Ma & Fry 
2000): even so, however, the previous reasoning would still 
be functional. The rmmcrical results for the vorticity. Sec. 4, 
also favor a ratio L/ro 0.1. 

This estimate is consistent with the fact that the dust 
model predictions agree with N-body simulations when the 
former are not UV-divergent, in some cases also up to the 
limit R^ ro (beyond the a priori expected range of validity, 
R ^ ro). as the example of the skewness in Sec. 3 shows, 
corrections are then of order {L/roY when pushing R oi ro, 
and thus possibly unobservable (e.g., the estimated absolute 
error of the measured skewness is about 0.1 (Juszkicwicz et 
al. 1993)). On the one hand, this good agreement of the 
convergent dust-model predictions supports the hypothesis 
of the mode-decoupling at the basis of the SSE. On the other 
hand, it suggests that even the dust-model UV-divergent 
perturbation results may provide good approximations when 
regularized in the framework of the SSE. 



6 CONCLUSIONS AND OUTLOOK 

I have explored the application of Eulerian perturbation the- 
ory to a recently proposed model of cosmological structure 
formation, the small-size expansion (SSE), which incorpo- 
rates corrections to the dust model. The SSE features a 
length scale L{t) which is identified with the typical size 
of the most recently formed virialized structures. The per- 
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turbation results suggest L/ro ~ 0.1 with ro the nonlinear 
scale. 

When the dust-model perturbation results are finite, 
the corrections are irrelevant. However, the length L{t) 
works as a natural UV-cutoff which regularizes perturba- 
tion expansions in the dust model which would be other- 
wise divergent. The SSE may also address multistreaming 
effects perturbationally, which are absent in the dust model. 
As an example, I have computed the vorticity generated by 
tidal torques: its relative importance is predicted to dimin- 
ish when the cosmic flow is observed on increasingly larger 
scales, and to be significant today on scales ~ 1 h~^Mpc 
with a value ^ Hubble constant. 

The SSE provides a new perspective and a possibly use- 
ful tool in the analytical understanding of the evolution by 
gravitational instability. The model based on the SSE must 
still be certainly subject to further investigation to check its 
internal consistency and eventual success: the results of this 
work can be contrasted to N-body simulations. As well as 
the predictions for the vorticity, it will be particularly inter- 
esting to compare regularized perturbation corrections, e.g., 
that for the density variance cr^. Finally, another interesting 
line of future research on the SSE is the application of the 
Lagrangian perturbation theory: this should help clarify the 
relationship to the truncated Zel'dovich approximation. 
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APPENDIX A: DERIVATION OF THE 
SMALL-SIZE EXPANSION (SSE) 

To make this work as self-contained as possible, I provide in 
this App. the derivation in Paper I leading to Eqs. (1). It is 
however more general in that a time-dependent smoothing 
length is allowed. 

The basic model is a system of nonrelativistic, identical 
point particles which (i) are assumed to interact with each 
other via gravity only, (ii) look homogeneously distributed 
on sufficiently large scales, which thus are assumed to evolve 
according to an expanding Fricdrnann-Lcmaitrc cosmologi- 
cal background, and (iii) deviations to homogeneity are rele- 
vant only on scales small enough that a Newtonian approx- 
imation is valid to follow their evolution. Let a{t) denote 
the expansion factor of the Friedmann-Lemaitre cosmolog- 
ical background, H{t) = d/a the associated Hubble func- 
tion, and Qb{t) the homogeneous (background) density on 
large scales, the comoving spatial coordinate of the 

a-th particle, Ua is its peculiar velocity, and m its mass. 
In terms of these variables the evolution is described by the 
following set of equations (Va denotes a partial derivative 
with respect to Xa)'- 

. _ 1 

a 



ila = Wa — HUa, 



(Al) 



Va • Wa = —4:nGa 

Vq X Wa = 0, 



— ^ S{Xa - Xff) - eb 



(A2) 



where Wa 'is the peculiar gravitational acceleration acting on 
the a-th particle. Finally, Eqs. (A2) must be subjected to 
periodic boundary conditions in order to yield a Newtonian 
description consistent with the Friedmann-Lemaitre solution 
at large scales (Ehlers & Buchert 1997). 

One can formally define a microscopic mass density field 
and a microscopic peculiar velocity field. 
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= (a:,t) 



a{tr 
{x,t) 



oc 



(A3) 



Coarsening cells of comoving size ~ L{t) are defined with the 
help of a smoothing window W{z) (see Paper I for details on 
the properties required on the window) . The coarse-grained 
fields associated to the microscopic fields are respectively 



g{x,t) = J^W 

Qu{x,t) = J ^ 



y\ 



w 



L 

\x -y 
L 



Qmiciy,t), 



(A4) 



Notice in particular that u is defined by a volume-average of 

the peculiar momentum density (i.e., by a mass-average of 
the velocity) and thus its physical moaning is transparent: 
u is the center-of-mass velocity of the coarsening cell. These 
fields obey the following exact evolution Eqs. (conservation 
of mass and momentum, following from Eqs. (Al)): 



dQ 

dt 



-3Hq- iv • (qu) + yV ■ 
d Li 



at a L 



(A5) 



New fields have arisen, whose physical meaning is also clear 
from its definitions: two vector fields, 



s(x,t) = j 
f{x,t) = J 



dy 

L3 

dy 



W 



W 



X 



■v\ 



L 

\x-y 



L3 ■■ \ L 
and two second-rank tensor fields. 



(y - x)Qmic(y,t), 

{y,t), 



mic ^rnic \ 



P{x,t) = j 

D{x,t) = 



Ay 

L3 

dy 

V- 



w 



w 



\x ■ 



■v\ 



L 

\x - y\ 



\irmc **rmc ""mzc 



Umic {y,t), (A6) 



(y - x)Q,nicU,nic{y, t) . 



(The definition of Wmic is Eq. (A3) with the replacement 
Ua Wa). If the evolution Eqs. for these fields are com- 
puted, new fields appear, and so on ad infinitum. A physical 
assumption is required to cut the hierarchy: the nonlinear 
coupling in Eqs. (A5) to the degrees of freedom below the 
smoothing scale L is assumed to be weak. To implement this 
idea mathematically, notice that Eqs. (A4) can be formally 
inverted to yield the following expansion in L (this is most 
Ccisily derived in Fourier space): 



Qmic — 



l__B(LV)'+o(LV)* 



l--B{LVy 



u, (A7) 



with the constant 



B = ^j dzz^ W{z) = ^ dz / W{z). 

These expansions just show how the microscopic fields can 
be recovered by taking into account higher and higher 
derivatives of the smoothed fields, i.e., finer and finer details 
in the spatial distribution. Now consider, e.g., the tensor 



field P, encompassing nonlinear mode- mode couplings: the 
weak-coupling assumption above then means that the dom- 
inant dynamical contribution of P to the evolution Eq. (A5) 
arises from the coupling between modes on scales > L. 
Hence, P can be replaced in this Eq. by the result of in- 
serting the expansions (A7) in the definition (A6): 

P ^ guu + BL'^g{d,u){d,u) + o{L V)^ 

Applying the same reasoning to the other fields: 

f ^ ew + BL^{Vq ■ V)w + o{L V)^ 

V -w = —4-wGa(g — Qb), 

V X w = 0, 



and 
s — »■ 



SLVe + o(LV)*, 



D ^ BL^V{gu)+o{LVf. 



(A8) 



Taking these results in Eqs. (A5), one obtains a formal ex- 
pansion in L (called the small-size expansion in the main 
text). Truncating after the second term, Eqs. (1) are recov- 
ered. 

The fields s and D have a purely geometrical origin 
and describe how the smoothed fields change simply be- 
cause the smoothing length is varied. Its diffusion-like ex- 
pression is thus barely surprising. In fact, truncating the 
expansions (A8) is exact for a Gaussian window, W{z) = 
exp(— 7ra^) and B = 1/2tv. In general, the approxima- 
tion (A8) is equivalent to replacing the arbitrary window 
W{z) by a Gaussian one with a width fixed by the constant 
B. 



APPENDIX B: VORTICITY VARIANCE 

In this App. I collect the mathematical details of the study 
of Eq. (18). 

We consider first the IR behavior. When P{k) ^ k" 
as fc — » 0, the integrals behave asymptotically in the limit 
g — » 0, {k/q) — »■ as (s := cos 0, := (fe, q)) 



|H(0,t)r dk k 



I' 

Jo 



l,"+2 



Jdq q"+^ J^ds s^(l-s'). 



which is finite and nonvanishing if n > —3, independently 
of the values of 7? and L. 

To address the UV behavior {k,q — »■ -|-oo), one makes 
the change of variable q ^ p = k — q. Assuming P{k) ~ fc" 
as k ^ +O0, we have asymptotically 



r 



H(fe^^)|^ / dppV-^^'" \W{Rp)fx 
Jo 



ds s^(l 



S ). 



The smoothing window over scale R renders the p-intcgral 
UV-convergent, no matter the value of L. From the def- 
inition of the function H, it follows that H(A;^,t) < 
S(0, t)e^^^*"'' (since L > by assumption), so that the 
fc-integral is also convergent as long as Lin 7^ 0, as expected. 
However, in the dust limit L = 0, the fc-integral converges 
only if n < -3/2. 
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In what follows, the particular case of an Einstein-de 
Sitter background is considered with the self-similar initial 
power spectrum 

dfe 



a two-dimensional integral in the finite domain < 2; < 1, 
< M < 2. 



P{k) 



-^n ' in. 1 



-X\w{k)Y 



(Bl) 



(27r)3 

where Tin = ro{tin). A Gaussian window is assumed, 
W{z) = exp(— TTZ^), and then 



In Eq. (18), three angular integrations can be performed 
immediately. Then, one exploits the symmetry of the inte- 
grand under exchange fe <-> q and changes to new variables 
k — » k/Lt, q — > kp: 



2Al 



2(n-|-3) 



(27r)4 ^2(n-H5) 



J dss^(l-s^^ 



/ dfefe^"+Md 

^0 Jo 



2x2 , 
-p ) X 



„ f k ps 



Inserting Eq. (Bl) in Eq. (11), one finds the linear dust 
prediction (i.e., up to corrections of order {L/Rf): ro{t) = 
rindl^''"'^^^ ■ As remarked in Sec. 4, a constant ratio A = 
Lt/ro{t) is taken. Then, the function H can be written in 
terms of the incomplete gamma function. However, it proves 
a better idea (particularly in order to compute au, numeri- 
cally) to perform first the fe- integral, and next, in the integral 
defining H, Eq. (19), to change variables t — » [ro{t)/rin]T = 



4/(n+3) 



2 r(n -1- 5) 



T. One finally gets: 
(5n + 23)at 



2"+3 



sr 



-| 2{n+3) 



1 + 



-(n+B) 



</' 

■'(ri, 



dri dr2 (T1T2) » H, 



with the function 



Tl -I- T2 

1 -h {RIU)' 



(B2) 



1 

dz z' 



n+2 



■ J ds s^(l 



^) 1 



-(n+5) 



obtained with another change of variable: z := 2p/(l +p'^). 
Only the domain < t/ < 2 is iutcrcsting. The function may 
exhibit a singularity as j/ — > 0"'" . A straightforward study of 
the convergence of the integrals near the limits s = 1, z = 1 
provides the asymptotic behavior 



1, 



y 



-3<n<--, 



< n. 



With this information, one can easily derive the scaling be- 
havior (20) from Eq. (B2). The long-time limit, t ^ tin, is 
obtained by sotting the lower limit of the r-intcgrals to zero 
(no divergence arises with the constraint —3 < n). In the 
exponent range —3 < n < 4, the long-time approximation is 
accurate to the 1% level for times such that at > 5. 

The integral (B2) is also well suited for numerical eval- 
uation. In fact, the s-intcgral in the function Hn can be 
computed analytically. And the introduction of the variable 
u := Ti+ T2 allows to compute one of the T-integrals analyt- 
ically too. In the end, the numerical computation reduces to 
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